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Abstract 

Combined with all density-dependent factors, the per capita growth rate of a species may be non- 
monotonic. One important consequence is that species may suffer from weak AUee effects or strong 
AUee effects. In this paper, we study the permanence of a discrete-time two-species-interaction model 
with non-monotonic per capita growth rates for the first time. By using the average Lyapunov functions 
and extending the ecological concept of the relative nonlinearity, we find a simple sufficient condition 
for guaranteeing the permanence of systems that can model complicated two-species interactions. The 
extended relative nonlinearity allows us to fully characterize the effects of nonlincarities in the per 
capita growth functions with non-monotonicity. These results are illustrated with specific two species 
competition and predator-prey models of generic forms with non-monotone per capita growth rates. 

Keywords: Permanence, Non-monotonic Per Capita Growth Rates, Alice Effects, Relative 
Nonlinearity, Two-Species Interaction Population Models 



1. Introduction 

This paper is concerned with the permanence of a discrete-time two-species interaction model with 
non-monotonic per capita growth rates that can be described by the following equations: 

xt+i = xtf{xt,yt) (1) 
yt+i = ytg{xt,yt) (2) 

where Xt and yt denote population densities of species x and y in season t respectively; f{x, y) and 
g{x,y) are per capita growth rates of these two species, which are nonnegative and twice differentiable 
in K^. In addition, at least one of f^i f^, ff ; f| change signs in R^, i.e., the per capita growth rates 
are non-monotonic. 

The coexistence of different species in ecological communities has been a main research theme in 
ecology. For deterministic models, the idea of permanence, which guarantees convergence on an interior 
attractor from any strictly positive initial conditions, is regarded as a strong form of coexistence. Perma- 
nence of dynamical systems has been studied by many researchers using Lyapunov exponents (Shreiber 
2000; Garay and Hofbauer 2003; Salceanu and Smith 2009) and average Lyapunov functions (Garay 
and Hofbauer 2003; Kon 2004; Kang and Chesson 2010) have used Lyapunov exponents and the notions 
of unsaturated invariant sets (i.e., a compact invariant set that is repelling respect to the system) and 
measures for Kolmogorov-type systems. Recent study by Kon (2004) uses an average Lyapunov function 
to show that nonexistence of saturated boundary fixed points is sufficient for permanence under certain 
convexity and concavity conditions on the per capita growth rates of the species. Kang and Chesson 
(2010) make use of the concept of the relative nonlinearity (Chesson 1994) to extend Kon's (2004) results 
beyond convexity and concavity conditions to arbitrary nonlincarities for two dimensional discrete-time 
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competition and prcy-prcdator models. However, Kang and Chesson (2010) as well as Kon (2004) make 
assumptions on the monotonicity of f{x,y) and g{x,y), e.g., ^ and ^ does not ehange the sign. In 
this article, we will drop this assumption and derive a easy-to-eiieek permanence criterion for a general 
two-species interaction model with the per capita growth rates being non-monotone by using the theory 
of average Lyapunov functions and extending the ecological concept of the relative nonlinearity. The 
original concept of relative nonlinearity is a species-coexistence mechanism that results from different 
species having different nonlinear responses to competition together with fluctuations in time or space 
in the intensity of competition (Chesson 1994& 2000; Kang and Chesson 2010). This extended concept 
allows us to fully characterize the effects of nonlinearities in the per capita growth functions, which are 
of major significance in the presence of fluctuated populations and non-monotonic per capita growth 
rates. 

The structure of the rest paper is organized as follows: In section 2, we introduce some basic termi- 
nology of semi-dynamical systems and prove the important lemmas that are critical to derive sufficient 
criterion for the permanence of two-species interaction models with non-monotonic per capita growth 
rates; In section 3, we give a simple sufficient condition for the permanence of prey-predation models 
with non-monotonic per capita growth rates; In section 5, we apply our results to specific competition 
models and prey-predator models with generic with non-monotone per capita growth rates. We conclude 
with a discussion of the broader implications and prospectus for future work in this area. 

2. Preliminary results 
Let 

X = {{x,y) : X > 0,y > 0}, 5, = {(.t,0) : a; > 0}, Sy = {(0,y) : y > 0} 

and 

S = S^[jSy, M ^X\S ^{{x,y):x>0,y>0}. 

Define a discrete two dimensional dynamical system H : X X where X is a metric space, and denote 
H^{x,y) = {x,y) and i/"(x-, y) = {xn,yn)- We call that H is positively invariant in X, if H^{X) C X 
for any n € Z_|_. We call that H is dissipative, if there exists _B > 0, for all initial condition {x,y) G X 
such that 

limsupmax{x„, j/„} < B. 

This implies that H has a compact absorbing set in X. For convenience, when a system H is dissipative, 
we will consider X is a compact metric space. 

We call that H is permanent in M, if there exists two positive numbers < 5 < -B, for all initial condition 
(x, y) G M such that 

b < liminf min{a;„, yn} < limsupmax{x„,y„} < B. 
This implies that H has a compact interior attractor B C M that attracts all the points in M. 

The boundary equilibrium {x*,y*) G S* is unsaturated if 

f{x*,y*)>landg{x*,y*)>l. 
The positive orbit of H with initial condition {xo,yo) € ^ is defined as 

7^(a;o,yo) = {{xi,yi) : {xi,yi) = W{xQ,ya), for all i G Z+}. 
The omega limit set of {x, y) G R\ is defined as 

Lo{x^ y) = {(^, 77) : H*^ (x, y) — > (^, rf) for some sequences tj — > 00}. 
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The omega limit set of a subset S G X of X is defined as 



The main goal of this paper is to find a sufficient conditions on / and g such that the system (l)-(2) 
is permanent in M when /(x, 0), /(O, y), g(a;, 0) and g(0,y) are non-monotonic in S, i.e., at least one of 
the follows change sign in 5", 

dfix,0) dfjO.y) dg{0,y) dgjx.O) 
dx ' dy ^ dy ' dx 

2.1. External Lyapunov exponents 

Let {(a;i,yi)}^Q to be the positive orbit with initial conditions {xQ^yo) e 5. Then the average per 
capita growth rates of species x,y with initial conditions (a;o,0) G Sx (or (0,yo) G Sy) after n — 1 
generations can be represented as 

.^(-0,0) = ffi^i^^ , .r(xo,0) = ffi^i^^ (3) 
rr(0,.o)=^S>^ , .-(0,yo) = S>Z(^ (4) 



oo 
n=l 



Define r'^'^(xo, 0), r^'^(a;o,0), r^^(0,yo), r^^(0, yo) as the hm sup of the sequences 

{r^ (2;o, 0)}- 1, {rr (^0, 0)}- 1, {rJl''(0, yo)}Jf=i, {r^^O, 2/o)} 

respectively. Moreover, we use 

f--(a;o,yo), r^"(xo,0), r^^a^cyo), ^^^0,2/0) 

instead if their limits actually exist. 

Notice that the quantity fy^{xo,Q) (or ra;y(0,j/o)) is the external Lyapunov exponent of S'^; (or S'y), 
which gives the average invasion speed of the invader y (or a;) (Rand, Wilson & McGlade 1994). If both 
f^^(a;o,0) and r^y(0,yo) are positive for all xq > 0,yo > 0, then species y and x are able to coexist. 
The interesting question is that what kind of conditions on f(x,y),g{x,y) can guarantee this, therefore 
guarantee the system (l)-(2) is permanent. In order to answer this question, we assume that f{x, y) and 
g{x,y) of (l)-(2) satisfy Condition H in the rest of the article, where H is stated as follows: 

H: Both f{x, y) and g{x, y) are strictly positive and twice differcntiable in X \ {(0, 0)} with /(O, 0) > 1 
and 5(0, 0) > 0. In addition, the following two limits exist and are continuous: 

lini /(x,y) = /(0,y) and linr 5(2;, y) = 5(2;, 0). 



Remark: Condition H implies that the system (l)-(2) is positively invariant in X. Now we need to 
show the following two lemmas first. 

Lemma 2.1. [Bounded population density] Assume that f{x,y) in the system (l)-(2) satisfies Con- 
dition H. If there exists < aoo < 1 such that 

lim sup f{x,y) = Oao, 

then the population density of species x is bounded by some positive constant. Moreover, if {xo, yo) & Sx 
with Xo > 0, then f^^{xo,0) = 0. 
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Proof. Define ai{x) — supygj{^ f{x,y), then the condition lim^j^^oo sup^gi^^ f{x,y) = Ooo < 1 indicates 
that for any e + a^o < 1, there exists a number large enough such that 

ai{x) < Ooo + e < 1 for all x > L^. 

Since ai{x) = supj^gg^ f{x,y), therefore, 

fix,y) < ai{x) < Qoo + e < 1 for all x > L^,y > 0. 

Now if an initial condition of species x is greater than Lg, i.e., xq > L^, then there exists a positive 
integer N such that xn < L^. Assume that this is not true, then for any positive integer n, we have 
Xn ^ L^. In particular, we have 

n-l 

Xn = Xn-if{xn-i,yn-i) = f{xi,yi) < a:o(aoo + e)" as n oo. 

This is a contradiction to the fact that Xn > La for all n e Z^. Therefore, there exists a positive integer 
N , such that xn < L^. 

Define Lm = -^e ma'X(2;,y)e[o,L5]2{/(a;, y)}. We claim that if xm < L^, then x„ < for all n > N. 
Suppose that this is not true, then there exists some positive integer P such that x^^p > Lm- Let 
Pm = min{p + : Xn+p > Lm}, then we have 

Xp^ ^ Pfn and Xp^ — \ ^ Lyy^. 

This implies that either 

Xp^-l < or < Xp^-i < Lm- 

If Xp^-i < Le, then 

Xprr, = Xp^-if{xp^-i,y) < if max {f{x,y)} = Lm 

(x,y)6[0,L,]2 

which is a contradiction to Xp^^ > Lm- 

If if < Xp^-i < Lm, then due to the fact that 

f{x, y) < Goo + e < 1 for aU (a:, y) e [L^, oo) x [0, oo), 

we have 

Xp^ = Xp^^if{xp^^i,y) < x-p,„_i(aoo + e) < Lm 
which is also a contradiction to Xp^ > Lm- Therefore, we have 

Xn < Lm for all n > N. 

This implies that for any initial condition (a;o, j/o) € X with a;o > 0, there exists a positive integer p,„, 
such that 

Xn < Lm for all n > Pm- 

Therefore, the population density of species x is bounded in the system (l)-(2). 

Next, notice that Sx is positively invariant, then for any initial condition in Sx, we have 2/„ ~ for 
all future n > 0, i.e., we have the following boundary dynamics, 

Xn+i = Xnf{xn, 0) for all n > 0. 

Then by applying Lemma B.l (Kang and Chesson 2010), we can conclude that for any initial condition 
{xq, 0) with xq > 0, the following inequalities hold 

< b < lim inf x'„ < lim sup x„ < Lm - 
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Let {{xi, 0)}^g to be the positive orbit 7+(xo, 0) starting at xq > 0, then we have 

. Inf . . In^ . In^ . In f 

— lim — — < lim inf < lim sup < lim sup — — — 

This impHes that for all > 0, we have 

r-(xo,0)= lim S>Z(^= lin, ^ = Q. 

Ti— >oo 77, n— >cxD 77 

Therefore, we have proved the statement. □ 

Remark: Lemma 2.1 gives an easy-to-check sufheient criterion for species x being bounded in a joint 
system (l)-(2), which can not only can apply to competition and prey-predator models but also can 
apply to the mutualism system. For example, we can apply Lemma 2.1 to show that the following 
mutualism model (5)-(6) (Kon 2004) is dissipative: 



xt+i = xte (5) 

i'22 

yt+1 = Vte (6) 



where < un < Vi2 and < V22 < i'2i- 

For convenience, define F{x, y) = In f{x, y), G{x, y) = In g{x, y) and Fi,Gi,i — x,y a,s the first partial 
derivative respect to i; Fn^Gu^i = x,y as the second partial derivative respect to i. Let {(xi,0)}^g to 
be a positive orbit 7"'"(a;o,0) with > 0. Then we have the following lemma: 

Lemma 2.2. [The external Lyapunov exponent] Let (x*,0) he a point in Sx such that f{x*,0) > 
and Fx{x*,0) 7^ 0. Then the following two cases hold if Condition H holds as well as 

lim f{x,0) = ai < 1. 

n— >cjo 

Case one: ///(0,0) > 1, then f^^{xo,0) = where Xq > 0. In addition, the external Lyapunov 
exponent of is 

fy^xo,0) = G{x*,0) - ^^F{x\0) + Ayixo). 

Case two: ///(O, 0) < 1, then f^^(a;o,0) < where Xq > 0. In addition, if p^'j^, < 0, then the 
external Lyapunov exponent of Sx satisfies follows 

f^^^ixo, 0) > Gix*,0) - ^M^F{x*,0) + Ayixo) 



whe 



Ay ~ lim sup 



dt 



and Xit ~ x* + (xi — x*)t. 

Proof. The condition that lim„_^oo fix,0) — ai < 1, indicates that the superior of the average growth 
rate of the species x without species y is less than or equal zero by applying Lemma 2.1, i.e., for all 
initial conditions xq > 0, we have 

~xx, r J2"=o ^Xn/Xo) ln(B/.To) 

r [xq, 0) — hmsup " ^ ~ limsup — ^ — < Imisup — ^ — ■ ~ 0. 
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Doing exact 2 order Taylor expansion on F{xi, 0), G{xi, 0) around x ~ x* gives: 



F(x„0) = F{x\Q) + F,{x\Q){x,~x*) + {x,~x* f ^l{l-t)F,.,{xu,Q)A 
G{x^,0) = G[x\Q) + G^{x\G){x, -x*) + [x^- x*f (1 - t)G^^{xu,0)dt 

where xa = x* + {xi — x*)t. Then we have 



<-(xo,0) 



This imphes that we have 



F,{x-\0) nF^(x*,0) 



This imphes that we can rewrite r^^(a;o) as foUows: 

ErJc/ G(3:, ,0) 



G{x\o) + Gx{x\o)^^^^^§^ + i:rJoH^.-^-)Vo^ji-*)G..(..,o)rft 



G(x*,0) + G,(x*,0) 



" rg°=(xo,0)-f(a;*,0) _ T,1:Zq (^i-^'f Sg il~t)F^^{xH,0)dt 
F^{x*,0) nF^ix'fi) 
T.7=oi^^-^'f lo {l-'t)G^U^.t,0)dt 
n 

r<(n,* C\\ J_ rl^{xo,0)G^{x' ,0) F(x-' .,0)G^{x' ,0) 
^(.^ ' ^> + F^ix'fi) F^ix',0) 

, ^-o\-^-'?!^ (i-*)[G..(^..,o)- ^-'-;:jr?ro';'"'" ] 



\dt 



Define Aj^(a:o) as 



Gx{x 



Fx{x*,0) 



Ay{xo) = (l-t) G^^{xit,0) 
Jo 

then we have 

Case one: If /(O, 0) > 1, then by applying Lemma 2.1, we have 

f^^(xo,0) = lim <-(xo,0) = 0. 

n— >c« 

This indicates that 



Case two: If < 0, then due to the fact that 



liminfr^^(a;o,0) < f^^(a;o,0) < 0, 



we have 



fy^ix,,0) > G(x%0)~ ^(-)°g-(f +liminf„^^ -"^-;-;)?-;-'-°) +A,(xo) 
> G{x*,O)~^i^i0^+Ay{xo) 

□ 

Remark: In the case that (a;*,0) is a nontrivial boundary fixed point, then F{x*,0) = 0, thus we have 
the following corollary from Lemma 2.2: 
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Corollary 2.1. [The external Lyapunov exponent] Assume that all conditions in Lemma 2.2 hold as 
well as F{x*,0) = 0, then 

Case one: ///(0,0) > 1, then the external Lyapunov exponent of Sx is 

r«"(a;o,0) =G(x*,0)+Ay(xo). 

Case two: If f {0,0) < 1 and '^p°\^, < 0, then the external Lyapunov exponent of Sx satisfies follows 

r^"GTo,0) >G(.T*,0) + A,(xo) 
where Aj, is defined as the same in Lemma 2.2. 

Remark: The expression Gxxix, 0) — '^"'"^p'^fl^^Q-l'^'^^ in Corollary 2.1 can be considered as the extended 
ecological concept of the relative nonlinearity from Kang and Chesson (2010). This allows us to give 
an easy-to-check sufficient criterion for permanence applicable to a broad range of situations and avoids 
checking the detailed information on u}{Sx) and calculating the external Lyapunov exponent f^^(xo,0). 
Depending on the signs of G{x*,0) and Ay{xo), there are four situations: 

Permanence: If both G{x*,0) and Aj,(a;o) are nonnegative and at least one of them is positive for 
all Xq e R, then the external Lyapunov exponent f^^(a;o,0) of Sx is positive. Thus we can apply 
Theorem 2.2 and its corollary 2.3 of Huston (1984) to show that species y is permanent. We will 
focus on this case in this article. 

Relative Permanence: Notice that it is possible that G{x*,0) < but f^^{xo,0) is still positive for 
almost every xq S Sx. This is the case when permanence fails due to the nontrivial boundary 
equilibrium point {x*,0) being saturated, which give a proper setting for the relative permanence 
(Kang 2011a; Kang and Smith 2011b). 

Boundary Attractor: The case when r^^(a;o,0) < for all xq G Sx represents the case when fluctu- 
ations associated with the non-point attractor undermine permanence because then the invasion 
rate is lower than that predicted by the point attractor (a;*,0). Thus, the system has no perma- 
nence due to existing attractors on the Sx- There are many models (Kon 2006; Kang et al 2008) 
presenting this scenario under some proper parameter ranges. 

Multiple Attractors: The case when f^^(xo,0) > for a dense set of xo € 5^; is a mixture of case 2 
and 3, which can generate rich dynamics such as riddled basin of attractions (Ashwin et al 1996; 
Ferriere and Gatto 1995; Kon 2006). 

3. Sufficient conditions for the permanence of a general two-species interaction models 

Let species x and y interact with each other in an ecology community and their population density 
can be described by the system (l)-(2). We arc interested in a two-specics system (l)-(2) that satisfies 
Condition H as well as the following conditions: 

Gl: There exists ai,a2 such that 

lim sup/(a;, y) = ai < 1 and lim supg(x,y) = a2 < 1. 

G2: All the nontrivial boundary fixed points are unsaturated, i.e., if (a;*,0) and (0,y*) are nontrivial 
boundary fixed points, then f{0,y*) > 1 and 17(0;*, 0) > 1. 
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G3: There exists some point (a;*,0) G such that ^;;(.t*,0) ^ and G{x*,0) - '^"^^^"if of > 
and the equaUty (7) holds 



Fy(,0,v*}G{Oy' 



G4: If .g(0,0) > 1, then there exists point (0,y*) G Sy such that Gy{Q,y*) ^ and F(0,y*) - 
> and the inequahty (8) holds 

'^yy^TV ) 

If g(0, 0) < 1, then < holds in addition to (8) holds. 

Condition H guarantees that the population of species x and y will not drop below and species x 
will not be too close to the origin (0,0). Condition Gl implies that both species x and y suffer from 
intra-competition, which drops their per capita fecundity below 1 if their population density is too 
large regardless of other species' population. The per capita growth rates /, g can be non-monotonic 
with respect to x, y. Condition G2 ensures that both species x and y have nonnegative growth rates 
at the nontrivial boundary equilibria. Condition G3-G4 guarantee that species x and species y have 
positive invading speed when their population density are rare. The last condition G4 also implies 
that species y is able to persist even if it is not persistent in its single state (i.e., g{Q,0) < 1). More 
importantly, (7)-(8) is an extended concept of the relative nonlinearity that introduced by Chesson 
(Chesson 2000). Relative nonlinearity is a species-coexistence mechanism that results from different 
species having different nonlinear responses to competition together with fluctuations in time or space in 
the intensity of competition (Chesson 1994). The expressions (7)- (8) can be treated as a general form of 
relative nonlinearity when Fx{x,0),Gx{x,0), Fy{0,y),Gy{0,y) are non-invertible (Chesson 2000). More 
generally, we can consider r^(a;) (or r^{y)) as a contribution to the invading speed of species y (or x) 
due to species x (or y) has fluctuated population in its single state, e.g., if species x (or y) has only 
point attractors (x*,0),i = 1..,m (or (0,y*),j = l,..,w), then r^(a;) = 0(or r^{y) = 0). If the species 
X has non-point attractors (i.e., fluctuated populations), then the contribution r^(a:) can be positive or 
negative. In this paper, we focus on the case when both r^(x) and r^{y) are nonnegative and the main 
goal of this section is to prove the following theorem: 

Theorem 3.1 (Sufficient conditions on permanence of two-species models). // the system (l)-(2) 
satisfies Condition G1-G4 as well as Condition H, then it is permanent in M. 

Proof. Since the system (l)-(2) satisfies Condition H-Gl, then according to Lemma 2.1, there exists a 
positive number L, such that for any initial condition {xo,yo) G X, we have 

limsupmax{x„,y„} < L. 

n— >-oo 

Therefore, the system is dissipative in X. In addition, according to Lemma 2.1, Condition H-Gl also 
implies f^^{xo,0) = for all xo > and the following two cases: 

1. fyy{0, yo) = for ah yo > if g(0, 0) > 1. 

2. fyy{0, yo) < for all yo > if 5(0, 0) < 1. 

Let {{xk,0)}'^Q and {(0, yk)}'kLo be positive orbits with initial conditions xo > and yo > respectively, 
and 

Xkt = X* + {xk - x*)t, ykt = y* + ivk - y*)t, fc = 1, • • • ,00. 
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Then apply Lemma 2.2, we have 

fy^xo,0) = G«,0) - + A,(xo) 

where 

Aj,(a;o) = Umsup 

and the following two cases depending on the value of g{0, 0), 

1. If g{Q, 0) > 1, then f-v{0, yo) = F(0, y*) - + A,(yo); 

2. If .g(0, 0) < 1 and < 0, then 



dt 



where 



^xiyo) = hmsup 



dt 



n 

Then according to Condition G2,G3 and G4, wc have 



1. If g(0,0) > 1, then 



inf f«^(xo,0) > and inf 7=^^(0, yo) > 0. 

xo>0 yo>0 



2. If g{Q, 0) < 1 and < 0, then 

inf f^''(xo,0) > and inf f^y(0,yo) > 0. 

a:o>0 yo>0 

Now we have the following two cases depending on the value of g(0, 0), 
1. If g{0, 0) > 1, then define P{x, y) = xy, thus 

p, s n-1 n— 1 

Therefore, wc have the following inequalities hold 

sup liminf > g(r-^^(.,o)+i„f.>o,v(.,o)) > J 
n>o{xom)eM^(x,o)es^ P[xQ,yo) 

sup liminf P(-'n,yn) ^ g(.,„(0,y)+i„f„>of.„(0,y)) >^ 

n>0 (^o,yo)eM^(0,y)fES„ P[xo,yo) 

In addition, for all (a;, y) g 5*, we have 

Pix,y) = 

Then by applying Theorem 2.2 of Hutson (1984) to the system (l)-(2), we can conclude that the 
system is permanent in M. 
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2. If ^(0,0) < 1, then define P{x,y) ^ x, thus 

n— 1 n— 1 



Pixo,yo) tn 



i=0 i=Q 

Therefore, we have the following inequalities hold 



sup 



lim inf 



n>0(xo,yo)&M^{y,0)£Sy P{xo,yo) 

In addition, for all {x,y) G Sx, we have 

Pix,y) = 

Then by applying Theorem 2.2 of Hutson (1984) to the system (l)-(2), we can conclude that species 
X is permanent. This implies that we can restrict the system in the space [b,Lm] x [0,Lm] with 
6 > 0. This allows us to apply Hutson's Theorem 2.2 (1984) again to obtain the permanence of 
species y in the joint system (l)-(2) by using the average Lyapunov function P{x,y) = y. Thus, 
the system is permanent in M . 

□ 

Remark: If the system (l)-(2) is discrete version of Loktcr-Volterra models, then F{x,Q)^G{x,Q) and 
F{0^y),G{0,y) are linear functions in x and y respectively. Therefore, Lokter-Volterra models are 
permanent if Condition G2 holds. 

Corollary 3.1 (Alternative version of Condition G3-G4). Assume that the system (l)-(2) satisfies 
Condition H, G1-G2 and Condition G'3-G'4, then the system is permanent in X, where Condition 
G'3-G'4 o.f^ stated as follows: 

G'3: There exists some point (x*,0) G Sx such that Fx{x* ^Q) ^ 0, G(x'*,0) — '^''^f ^ ^ '^'^'^ 



for all Xi > 0, we have 



ry{x,)^ / (l-t) 



Gxxi^it, 0) 



Gxix*,0)Fxxix,t,0) 



dt > 0, where xa — x* + {xi — x*)t. (9) 



Fx{x*,0) 

G'4: // .9(0,0) > 1, then there exists point {0,y*) € Sy such that Gy{0,y*) ^ 0, J^(0, y*)- 



and the inequality (10) holds for all yi > 0, 



> 



r^y,)^ / (l-t) 
Jo 

//.g(0,0) < 1, then 



Fyy{0,yit) 



Fy{Q,y*)Gyy{Q,yu) 



G,(o,r) 



dt > 0, where yu = y* + (yi - y*)t. (10) 



^^("'^*^0a.dF(0,,-)-^^^»:^:)^(.^ 



G,(o,r) " ' G,(o,r) 

for some point (0,y*) £ Sy in addition to (10) holds. 

Remark: The proof of Corollary 3.1 is similar as the proof of Theorem 3.1 and is straightforward, 
therefore we omit the details. In the case that (a;*,0) and (0,2/*) are nontrivial boundary equilibria 
points (i.e., F{x*,Q) = G{0,y*) = 0), then the conditions F{0,y*) - ^''^^^(o^l"''^'^ > and G{x*,0) 
^-^^;gfg''°) > can be dropped. 
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4. Sufficient conditions for the permanence of a general prey-predator model 

Let species a; be a prey and y be the predator of x in an ecology community and their population 
densities can be described by the discrete system (l)-(2). In a prey-predator model, the predator per 
capita growth rate is an increasing function of prey abundance and its per capita growth rate is commonly 
negative at zero prey abundance. Thus, the predator cannot survive without the prey, i.e. g(0,?/) < 1 
for all y > 0. We are interested in a prey-predator model (l)-(2) satisfies Condition H as well as the 
following conditions: 

PI: > and g{0, y) > for aU (x, y) E X. 

P2: There exists ai,a2 such that 

and 



lim sup/(x,y) = oi < 1 



for any a; > 0, lim g{x, y) = 02 < 1. 

P3: All the nontrivial boundary fixed points are unsaturated, i.e., if (x*,0) is an nontrivial boundary 
fixed point, then g{x*,0) > 1. 

P4: There exists some boundary fixed point (x*,0) such that 

rvix) = G,A^.,0)~^^^l^^^^^^p^>0,fora\\x>0 (11) 
or there exists some nontrivial boundary fixed point (a:*,0) such that for all Xi > 0, we have 



ry{x) = (l-t) 



Gx{x 







Gxxixit,0) 



F4x*,0) 



dt > 0, where xu = x* + {xi — x*)t. (12) 



Condition P1-P2 and H guarantees that the prey x is permanent without the predator y. InPl, the 
condition g{0,y) > y > implies that the predator has positive per capita growth rate; In P2, the 
condition limy_i.oo gix, j/) = 02 < 1 for any given a; > implies that for any given prey population size, x, 
the finite rate of increase of the predator drops below 1. Such behavior is normally described as predator 
interference, i.e. predator individuals interact negatively with one another limiting their ability to hunt 
prey. Condition P3-P4 makes sure that predator y is persistent with presence of prey x. We want to 
show the following theorem: 

Theorem 4.1. If the system (l)-(2) satisfies Condition H, PI, P2, P3, P4, then it is permanent 
in M. 

Proof. Since the system (l)-(2) satisfies Condition P1-P2, then by applying Lemma 2.1, we know that 
there exists a positive number L„i, such that for any initial condition {xq, y^) X , wc have 

limsupx,i < Lm and f^^{xo, 0) = for aU xq > 0. 

n— fcxD 

This implies that for any {xo,yo) G X, there exists a positive integer iV, such that 

Xn ^ Ljn for all n > N. 

By using the fact that g{x,y) is increasing with respect to x, we can conclude that there exists N large 
enough such that 

Vn+l = VngiXn, Vn) < 2/nff(im, Vn), for all n> N. 
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Since and for any given x > 0, we have 



lim g{x,y) = 02 < 1. 



Therefore, we can apply Lemma B.l (Kang and Chesson 2010) to conclude that there exists a positive 
number B > such that 

lim sup y„ < B. 

n— ^00 

Therefore, the system is dissipativc in M. 

Now we can show that prey x is persistent in the joint system (l)-(2) by applying Hutson's Theorem 
2.2 and its Corollary 2.3 (1984) with an average Lyapunov function P{x,y) = x. Then, we can restrict 
the system in the space [b,L,n] x [0,-B]. This allows us to apply Hutson's Theorem 2.2 (1984) again 
to obtain permanence of y in the joint system by using the average Lyapunov function P{x, y) = y. 
Therefore, the system is permanent in M . Hence, we have shown the statement. □ 

Remark: Note that ^^'^q^^ > 0, thus for any nontrivial boundary fixed point (x* , 0) , we have Gxix* ,0) > 
0, thus there are the following two simplified cases: 

1. If both F{x,0) and G{0,x) are convex (or linear), then Fx{x*,0) < and G{x*,0) > indicates 
that the system (l)-(2) is permanent. 

2. If i^(x, 0) is convex (or linear) and G{x, 0) is concave (or linear), then Fx{x* ,0) > and G{x* , 0) > 
indicates that the system (l)-(2) is permanent. 

5. Applications 

In this section, we apply the results to particular competition and prey-predator models. 

5.1. A competition model with strong Allee effects in one species 

Multiple stable states occur when more than one type of community can stably persist in a single en- 
vironmental regime. Simple theoretical analyses predict multiple stable states for single species dynamics 
via strong Allee effects. Perhaps the most common Allee effect occurs in species subject to predation by 
a generalist predator with a saturating functional response (Schreiber 2003). The population dynamics 
of a species suffering Allee effects due to predator saturation can be modeled as 

Nt+i = Ate'-(i-^'/^)/(7Vt) (13) 

where I{N) = e~ represents the probability of escaping predation by a predator with a saturating 
functional response and m, s represents predation intensity and the proportion to the handling time 
(Hassell et al. 1976) respectively. Non-dimensionalizing (13) by setting yt = Nt/K and h — sK gives 

yt+i = yte'-^'-y^^-^ (14) 

whose dynamics depends exclusively on the quantities r; b; and m. If r < m < ''^''^"'"^ and 6 > 1, 
then (14) suffers from strong Allee effects which leads to two positive equilibria that can be expressed 
as follows: 

r{b - 1) ± Jr'^il -I- 6)2 - Ambr 

V = TT ■ 

2rb 

Assume that the population dynamics of species x and y can be modeled as follows: 

xt+i = xte'^'^^-^'^-^^^ (15) 



12 



where r, , i ^ 1,2 represents the intrinsic growth rates of species x and y respectively; the term e~ i+^'n't 
represents that species x suffers the saturated inter-competition from species y; the term e-aa;(a;-c) 
represents that species y benefits from species x if the abundance of x below the threshold c, otherwise, 
species y suffers inter-competition from x. This system (15)-(16) models the following two features of 
the interactions between two species: 

1. Multiple equilibria: the species y suffers strong Allee effects in the absence of species x, i.e., 
^2 < ^2 < and 62 > 1. 

2. Threshold: If the population of species x is below c, then species y benefits from the presence 
of species x\ however, if the population of species x is above c, then species y suffers from inter- 
competition from X. 

It is easy to check that the system (15)-(16) satisfies Condition H-Gl. The boundary fixed points of 
the system are 

(0,0),(l,0),(0,yi*) and (O,^;) 

where 

, r2(62 - 1) - v/^f(TT6^?~"4^^^^A^ , , ^2(62 - 1) + yjrl{\ + h^Y - ^m^h^r-i 

Vi = ; and 2/2 — ; • 

2r202 2r202 

Thus, if r-i — + a(c — 1) > and r\ — > 0, then the system (15)-(16) satisfies Condition G2. 

Let 

F{x,y) = ri{l - x) - and G{x,y) ^ r2{l - y) - ——: ax{x-c). 

1 + biy 1 + b2y 

Then 

2bjjrii 



F,(x,0) = -ri, F,,(x-,0)-0, F,(0, y) = - j^:^, FyyiO,y) 



(l + fcia)3 

G,(x,0) = a(c~2x), GrA^,0) = ^2a, Gy{0,y) = -r^ + Gyy{Q,y) ^ 

Notice that 

F(0,0)=ln/(0,0)>0, G(0,0)=lng(0,0)<0, ^Ml = <o 

Gy(0,0) 02»Ti2 - r-2 



and 



j,^QQ^ F,(0,0)G(0,0) mi(m2-r2) 



Gj;(0,0) 52m2-r2 ' 

then the system satisfies Condition G4 (or G'4) if ri > mi and 

Fy{Q,Q)Gyy{Q,y) _ 2b^m^ 2mi6|m^_ 

'"'^ 'y' G,(0,0) (l + 6iy)3 (fe2m2-r2)(l-l-52y)3 

This implies that if 

hi < 62 and 7711(62^2 - > b\'m2, 
then the system satisfies Condition G4. In addition, we have 



r^(xO = / (1-t) 



Gxx{xit, 0) 



G,(l,0)i^,,(a;,t,0) 



dt 



Since < .t,; < ^-^q — , i = 1, • • • , 00, thus if 

G(l, 0) + (x, - lfry{xi) > a{c - 1) + r2 - m2 - max{a, a i— 1 j } 

holds, then the system satisfies Condition G'3. Therefore, based on the discussion above, we can apply 
Theorem 3.1 and Corollary 3.1 to gain the following corollary. 
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Corollary 5.1. The system (15)-(16) is permanent in M , if the following conditions hold 
Condition 1: r2 < m2 < ^^^4^^^^ o-JT-d 62 > 1 

Condition 2: a{c — 1) + r2 — TO2 ~ max{a, a (^^—f^ — - } > '"^'^ ''1 > max{TOi, ^} 

Condition 3: bi < 62 cind mi(&2n^2 — '"'2) > b\m2 
Remark: Examples of parameter's values that satisfy conditions in Corollary 5.1 are 

a = ?7i2 = 5i = 1, 62 = 2.5 > 1, TOi = 4, r2 = 0.85 < 1, rl = 4.1, c> 21. 

5.2. A competition model with weak Alice effects in both species 

A species suffering weak Allee effects due to predator saturation can be modeled by (14) when r > m. 
Then a two species competition model that is subject to weak Allee effects due to predator saturation 
for both species can be modeled as 

xt+i = xte''^^'-"*^"^^-'^!?/* (17) 
yt+i = yte^^^'-y*^-T^-^^^^ (18) 

where a^, i ~ 1,2 arc parameters that measure inter-competition between two species. If > mi,i — 1, 2, 
then the boundary fixed points of the system (17)- (18) are 

(0,0),(x*,0) and (0,y*) 

where 



, _ ri(&i - 1) + y^rfil + biY - Amibin , _ ^2(62 - 1) + ^ rl{l + 62)^ - 4771262^2 

X — Rncl y - — - . 

2riOi 2r202 

It is easy to check that the system (15)-(16) satisfies Condition H-Gl. Let 

F(x,y) = ri(l - x) - aiy and G{x,y) = r2(l - y) - 02- 

1 + bix 1 + 622/ 

Then F(0,0) = ln/(0,0) > 0, G(0, 0) = ln.g(0,0) > and 

F,{x,0)^-n + j^^, F,4x,0) = -j^^, Fy{0,y) = -a^, Fyy{0,y) = 

G4x,0) = -a2, G.4x,0) = 0, G,(0, y) = -r2 + j^^, GyyiO,y) ^ -j^^ 



Notice that 



and 



Gxx{x, 0) 



Fy{Q,y*)Gyy{Q,y) 2a,blm2 



G,(a;*,0)F,,(x,0) _ 2a2fo?TOi 



thus, if r2 > (i^yjy.)2 and n > jf^^i^^, then the system (17)-(18) satisfies Condition G3-G4. In 
addition, if F{0,y*) > and G{x*,0) > 0, then the system (15)-(18) satisfies Condition G2. Therefore, 
based on the discussion above, we can apply Theorem 3.1 to gain the following corollary. 

Corollary 5.2. The system (17)-(18) is permanent in M , if the following conditions hold 
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Condition 1: Vi > mi,i ^ 1,2 

Condition 2: F{0,y*) > and G{x*,0) > 

Condition 3: r-i > ff"",,, and n > nlT\,2 

Remark: Notice that F^ix, 0) and Gy{0, y) will change signs as x and y increasing to certain thresholds. 
Thus, permanence theorems in Kang and Chesson (2010) and in Kon (2004) fail while Theorem 3.1 
applies. 

5.3. A prey-predator model 

One-dimensional heuristic models of single-species population dynamics subject to the weak AUee 
effect can be represented as (19) (Lewis and Kareiva 1993; Amarasekare 1998a&1998b; Keitt et al. 2001) 

^ = x{r^h{x~af) (19) 

where all the parameters are strictly positive and r > ba^ . Then the discrete version of (19) can be 
represented as 

xt+i = Xte('-^("'-")') (20) 

Let Xt and yt represent the population density of prey x and predator y at generation t respectively. 
Then a prey-predator model with prey subject to weak Alice effects can be defined as 

xt+i = a;te(''-^(^'-")')-"i^' (21) 

yt+i = yte-'-^-^y' (22) 

where all the parameters arc strictly positive and r > ha} . The fact that all the parameters are nonneg- 
ative implies that this prey-predator system (21)-(22) satisfies Condition H, PI, P2. Thus, (21)-(22) 
is dissipative. Since r > ba^, then the only nontrivial boundary equilibrium is 

(a;*,0) = {a + y^,Q). 

Therefore, the system satisfies Condition P3 since 

5(a;*,0) = e"^("*)' > 1. 

Let 

F{x, y) = r — b[x — a)^ — ciy and G{x, y) = C2X^ — dy. 
Then F(0, 0) = In /(O, Q) = r-ba^ > 0, 0(0, y) = In 5(0, y) ^-dy<0 and 

F4x,0) = -2b{x-a), F,4x,0)^-2b, Fy{0,y) = -ci, Fyy{0,y)^0 
G^{x,0)^2c2x, G,,(x,0) = 2c2, Gy(0,y)==-d, GyyiO,y)^0 

This implies that the system satisfies Condition P4 if 

rvix) = GU-, 0) - ^j^FU-, 0) = 2c. - = 2c,(l - ^±^) > 

Fx{x*,0) b{a-x*) y/rb 

Therefore, we have the following corollary by applying Theorem 4.1: 

Corollary 5.3. Assume that all the parameters of the system (21)-(22) are strictly positive and 
r > ba^ . Then the system is permanent in M , if the following conditions hold 

1- °+V^ >o ( i.e.^a' <r{b + l/b-i)). 



Remark: Notice that F^(x,0) will change signs as x depending on the value of x. Thus, permanence 
theorems in Kang and Chesson (2010) and in Kon (2004) fail while Theorem 4.1 applies. 
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6. Discussion 



Population per capita growtli rate describes tire per capita rate of growtli of a population, as the fac- 
tor by which population size increases per year, conventionally given the symbol A — —^f^, or as r = In A. 
Therefore, population per capita growth rate is the summary parameter of trends in population density 
or abundance, which tells us whether density and abundance are increasing, stable or decreasing, and 
how fast they are changing (Sibly and Hone 2002). The per capita growth rate of a population can be 
broken down into negative density-dependent, density- independent, and positive density-dependent fac- 
tors (Shreiber 2003). Negative density-dependent factors include resource depletion due to competition 
(Tilman 1982), environment modification (Jones et al. 1997), mutual interference (Arditi and Akcakaya 
1990) and cannibalism (Fox 1975). Positive density-dependent factors include predator saturation, co- 
operative predation or resource defense, increased availability of mates, and conspecific enhancement of 
reproduction (Courchamp et al.l999; Stephens and Sutherland 1999; Stephens et al. 1999). Combined 
with all density-dependent factors, population per capita growth rate can be negative density-dependent 
in some range of population density and positive density-dependent in some other ranges, i.e., the per 
capita growth rate of a population may be non-montonic with respect to its population density. 

Processing multiple nontrivial equilibria is one important consequence of a species with non-monotonic 
per capita growth rate. The possibility that plant and animal populations have multiple positive equilib- 
ria has received considerable attention in the ecological literature. Theory and observation indicate that 
natural multi-species assemblies of plants and animals are likely to possess several different equilibrium 
points (May 1977). Ecological examples include fish (e.g., Peterman 1977; Spencer and Collie 1997), 
insects (e.g., Ludwig et al. 1978; Kuussaari et al. 1998; Solow et al. 2003), and phytoplankton (e.g., 
Beltrami 1989). Moreover, subtidal marine ecosystems in general, and reefs in particular, have several 
attributes which favor the existence of multiple stable states (Knowlton 1992). Beyond its scientific 
interest, this possibility also has important implications for the conservation and management of natural 
systems (Carpenter 2001). In this article, we focus on the permanence of a general two-species inter- 
action model with non-monotonic per capita growth rate for the first time. Our results broaden the 
applications of permanence for general two-species interaction models. 

Many mathematicians (Fonda, 1988; Freedman & So 1989; Shreiber 2000; Garay and Hofbauer 2003; 
Kon 2004; Salceanu and Smith 2009; Kang and Chesson 2010), have studied sufficient conditions on 
permanence of the dynamical systems. In this paper, we give an easy-to-check sufficient condition for 
guaranteeing the permanence of systems that can model complicated two-species interactions by using 
the average Lyapunov functions and extending the ecological concept of the relative nonlinearity. The 
extended relative nonlinearity allows us to fully characterize the effects of nonlinearities in the per 
capita growth functions with non-monotonicity. These results are illustrated with specific two-species 
competition and predator-prey models. In particular, we would like to point out that our result can 
apply to these models while theorems in Kang and Chession (2010) as well as Kon (2004) fails. In 
addition, our theorem can be easily extended to the mutualism models like (5)-(6). Extend our results 
to higher dimensional discrete-time population models (i.e., food- web dynamics in ecology) will be our 
future work. 
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